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Topology is familiar mostly from mathematics, but also natural sciences have found its concepts useful. 
Those concepts have been used to explain several natural phenomena in biology and physics, and they are 
particularly relevant for the electronic structure description of topological insulators and graphene systems. 
Here, we introduce topologically distinct graphene forms - graphene spirals - and employ density- functional 
theory to investigate their geometric and electronic properties. We found that the spiral topology gives rise 
to an intrinsic Rashba spin- orbit splitting. Through a Hamiltonian constrained by space curvature, 
graphene spirals have topologically protected states due to time-reversal symmetry. In addition, we argue 
that the synthesis of such graphene spirals is feasible and can be achieved through advanced bottom-up 
experimental routes that we indicate in this work. 
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In mathematics, topology analyzes how the properties of objects preserve under continuous deformations. But 
the interest to topological analysis is not restricted to mathematics alone; it spans also through biology, 
chemistry and materials science. In protein systems topology determines when protein folding sustains the 
rest of their cellular life\ In condensed matter physics, topology dominates several quantum phenomena, such as 
quantum-HalP, spin-HalP, and Aharonov-Bohm^ effects, as well as the physics of topological insulators^ 

In topological insulators, the surface electronic states are governed by topological features, making their 
quantum information robust against impurity scattering. Such robustness, by being protected by time- 
reversal-invariant Hamiltonian, could pave a reliable avenue toward fault- tolerant quantum-computing techno- 
logy^. Experiments via angle- resolved photoemission spectroscopy performed in Bi2Se3 compounds^ and 
Bii_;cSb;c alloys^ have shown signatures specific to topological insulators, such as large bulk energy gap and a 
single-surface Dirac cone associated to its topologically protected state. Dirac cones make the physics of graphene 
and topological insulators similar, even though graphene has two Dirac valleys with spin degeneracy while 
topological insulators have only one Dirac valley without spin degeneracy^. In addition, graphene can exhibit 
topologically protected quantum-Hall states with applied perpendicular and periodic magnetic fields^°. 

In the absence of structure inversion symmetry, surface states may split because of Rashba spin-orbit inter- 
action^ \ This splitting has been verified in thin films and semiconductor heterostructures having an inversion 
asymmetry of the confining potentials^^"^^. The splitting arises from the Rashba interaction Hamiltonian^ 
HR = a[Exp) •?, where a is the Rashba coefficient, s is the spin of an electron moving with momentum p = fik 
in an electric field E^^. In a two-dimensional non-interacting electron gas, therefore, Rashba spin-orbit interaction 
splits the parabolic energy bands in two, e+ (k) =Pt^k^ /2m* ± a/c, where is the effective mass^^. Even though 
spin-orbit interactions are often intrinsic, such as the spin-orbit-induced —0.1 meV energy gap in graphene^^, 
Rashba splitting is interesting for applications because the control over an external electric field makes it extrinsic 
and tunable^^. 

Extrinsic spin- orbit manipulation in graphene has been probed by external electric fields^°, by doping^ ^ by 
mechanical folding^^ and by depositing graphene on substrates^^. Especially substrate interfaces, by always 
involving inversion asymmetries, strengthen Rashba interaction^^. In graphene-based structures, however, the 
origin of the Rashba interaction is unlike in any other material: the interaction, although similar to the usual spin- 
orbit effects, arises not from the real electron spin, but from the spin related to the two non- equivalent atomic sites 
in the unit cell, the pseudospin. In graphene the Rashba splitting occurs around /c-points displaying time-reversal 
symmetry^^, meaning that the splitting is seen around Brillouin zone centers and zone boundaries'^. In graphene 
the zone boundary points K and K' are non- equivalent, but they have the same energy and they are connected by 
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time-reversal symmetry^^"^^. This symmetry makes the connection to 
topological insulators. The symmetry protects a pair of gapless hel- 
ical edge modes in topological insulators belonging to Z2 class^^. In a 
prototype of such topological insulator, an intrinsic spinorbit coup- 
Ung induces a topological mass term in the electronic structure of an 
atomic hexagonal frame^^'^^. Also, it is important to mention that 
most often Z2 class materials are discussed, at least experimentally, 
including heavy elements with strong spin-orbital coupling like 
Bi2(Se,Te)3 compounds^. 

Till now spin- orbit effects with highly protected topological edge- 
states in graphene have been obtained only by external perturbations 
such as external fields, heavy mechanical distortions, or chemical 
doping. In this letter, therefore, we investigate the above-discussed 
topological signatures in distinct graphene systems: the graphene 
spirals (see Fig. 1). As a central result, we find out that the spiral 
topology creates an intrinsic Rashba splitting. While in canonical 
illustrations of the effect in which an external magnetic field is 
applied to induce electron precession, in graphene spirals the track 
in which electrons move already displays an helical topology; elec- 
trons are constrained to move along an helical path in k- space and, 
due to unbreakable structure inversion symmetry, the Rashba-like 
band structure topology becomes an intrinsic effect for this material 
class. Our results demonstrate that graphene spirals naturally possess 
robust topological states as those observed in topological insulators. 

Results 

Graphene spirals are distinct from the helical graphene motifs 
reported earUer, such as graphene stripes and ribbons bent to 
spiral-like shapes^^"^^. In those earlier motifs, the starting point has 
been a graphene ribbon itself, with regular edges and curvature- or 
strain-modified 7i-electron system^^'^^. In such systems there has 
always been one-to-one mapping between the helical structure and 
flat graphene. In other words, it is always possible to build such 
systems by cutting them out of an infinite and flat graphene sheet. 
On the contrary, in graphene spirals such a mapping does not exist 
(see Fig. 1). Spirals are one- dimensional systems, while they still have 
a graphite layered structure containing perfect hexagons. Edge pro- 
flles alternate between armchair and zigzag shapes. Since spirals' 
local structures resemble graphene, they facilitate chiral topo- 
logy without overly perturbing the 71 -electron system. The largest 




(c) (C) , 




Figure 1 | Examples of graphene spirals, (a) Right-handed spiral, 

(b) interconnected right- and left-handed double spiral, and 

(c) loxodrome-like spiral, (a-c) Mapping of the (a-c) structures into helical 
stairs to highlight the topology of their curved space. 



perturbations take place at the inner edge of spirals where the strain 
is largest. 

First, we illustrate a notation scheme that we elaborated to identify 
the graphene spiral models (see Fig. 2). This notation is based on the 
number of armchair (ac) and zigzag (zz) segments used to complete 
a full coil of the spiral along the inner and outer edges. We estab- 
lished that the spiral axis goes through the center of a hexagon where 
7i/6- symmetry lines merge. Distinct spiral classes can be formed by 
combining ac and zz segments. The logic of our notation scheme can 
be better illustrated by projecting the spiral systems onto the plane. 
According to our construction rules, the spirals can only be formed 
by hexagons, their edges can simply contain ac- and zz-fragments 
and there must be no untangled bonds. How wide is the spiral can be 
determined by the outlined circles with radii rj, r2, r/ being Yj = 
1.23 + 2.46; A. These circles always enclose the flrst series of ac 
segments crossing the symmetry line SI. The remaining of the spiral 
is completed by as much zz segments required to maintain the 
perfect hexagonal frame. In this way, we define the notation 
diC[m]zz[nin-nout]y where m is the number of ac units and 

^in/out IS 

the number of zz units used to complete, respectively, the spiral inner 
and outer edges. Spirals with ac segments which are flipped towards 
the origin have an additional in their notation, e.g. ac*[2]zz[l-3] 
[see Fig. 2 (c) and (d)]. These conventions impose the spirals to be 
invariant under an axial translation of b (interlayer separation), 
under an axial translation of b/6 combined with a rotation of 7i/3. 
Furthermore, it is worth mentioning that almost any parametric 
curve can be used to represent edges formed under such hexagonal 
basis, however analysis involving more complex edge topologies 
are far beyond the scope of this work. The geometry and electronic 
structure of the built graphene spirals were investigated within 
density functional theory (DFT) implemented within SIESTA 
package^^'^". To confirm our results, SIESTA calculations were com- 
pared to other methods such as single Ti-band tight-binding (see 
Supplementary material), density-functional tight-binding imple- 
mented within DFTB+ code^^'^^, and VASP density-functional pack- 
age'^l Detailed description of the used parameters and calculation 




Figure 2 | Notation of graphene spirals edges: Right, Middle and left 
panels depict three-dimensional (3D) top view, conventional, and 
mirrored spirals with flipped ac segments, respectively. Note that on the 
right panels only carbon-carbon bonds along the inner and outer edges of 
the spirals are highlighted. Armchair (ac) segments are marked in red while 
zigzag (zz) segments are drawn in green colour. 81, 82 and 83 black lines are 
the hexagonal symmetry axes (71/6 wedge angles) and pass through the 
origin point located on the center of the spiral. The blue circles outline 
possible extensions that the spirals can hold along their inner and outer 
edges, (a) 3D view of ac[ l]zz[0-4], with its planar projection shown in (b). 

(c) Planar projection of ac* [ l]zz[0-4] spiral where ac segments are flipped. 

(d) 3D view of ac*[2]zz[l-3], with its planar projection shown in (f). 

(e) Planar projection of ac[2]zz[l-3] spiral. 
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conditions are presented in the end of the manuscript. We obtained 
the band structures for distinct graphene spiral models which had 
their atomic configuration fully optimized. 

Ab initio band structures obtained for graphene spirals of different 
widths are shown in Fig. 3. One can see that the two narrower spirals 
are semiconductors with an energy gap of approximately 2.3 and 
1.0 eV for ac[l]zz[0] and ac[l]zz[0-l], respectively. Differently, 
ac[l]zz[0-2] and ac[l]zz[0-3] reveal semimetallic states. However, 
more intriguing than the change of electronic character in response 
to geometrical variations is the peculiar band splitting at the Brillouin 
zone center. These band splittings represent the Rashba effect, the 
topological signatures in the electronic structure of graphene spirals. 
This is our central result. The splitting is robust and present in all 
spirals. The origin for the Rashba effect is related to how the Dirac 
particles in graphene couple with the spiral curved space. Graphite 
crystals, for instance, already manifest a natural splitting of the n- 
bands close to the Fermi level due to the two non- equivalent carbon 
atoms in the unit cell. Since the number of atoms in the spiral unit 
cell - assumed to be an unique coil - is considerably higher, one 
expects a superior number of splittings in these systems. In addition, 
the interlayer interaction in spirals also plays an important role. The 
usual graphite layer stacking provides an uniform potential profile 
over the whole sample which is not the case for the spirals. They 
possess chiral symmetry which leads to electron-hole symmetry 
breaking. Under chiral symmetry operation SJ, the wavefunction 
composed of a set of molecular orbitals within the unit cell get trans- 
lated (by t) and rotated (by 0) simultaneously^^. The operator Sq 
replaces the translation operator T in Bloch's theorem^^ It is there- 
fore important to distinguish reciprocal k-vectors of the linear (k) 
and curved or chiral (k') systems. These vectors are related by k' = 
{hlz)k, as can be demonstrated by applying S consecutively until a full 
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Figure 3 | Band structures of graphene spirals, (a) ac[l]zz[0], 
(b) ac[l]zz[0-l], (c) ac[l]zz[0-2], and (d) ac[l]zz[0-3]. Fermi energy is set 
at 0 eV. Vertical axes span the /c- vector intervals [-n/byii/b] . Narrow spirals, 
semiconductors with 2.3 - 1.0 eV energy gaps; wide spirals, semimetals. 
Bottom panels depict a zoom over low energy bands around the P. The 
band closest to the Fermi energy is highlighted in red while the second 
closest is colored in blue. 



turn is completed. The enlarged Brillouin zone in chiral systems is 
often referred to the Jones zone QZy^. It can be shown that in a 
constantly curved space systems, a bisector reduction of whole the 
Brillouin space must have two components as for any system 
described by a Hamiltonian under chiral symmetry. 

Let us now illustrate in detail the origin of the Rashba effect in our 
graphene spirals by considering a linear chain model as starting point 
(see also Supplementary material). The chain periodicity is given by 
the unit cell length a, and chain's one-dimensional potential obeys 
V(z) = V(z + a). When the unit cell gets more atoms, the reciprocal 
space shrinks, the bands fold and their number increases. For unit 
cells with an even number of atoms, two types of bands exist: bands 
crossing the Fermi level (i) at the T point and (ii) at the edges of the 
Brillouin zone. At the T point, the energy states are quantized as for a 
finite ring with M atoms, 6n{k) = — It cos{2nn/M), where n = 1,. . ., 
M and t is the hopping parameter. Only when 2n/M = ±l/2,±3/ 
2,±5/2,. . ., a band will cross the Fermi level at T point. For unit cells 
with an odd number of atoms, the band crossing is shifted from T 
point and an energy gap opens. A reminiscence of this cosine trend 
exists also in polyacetylene-like spiral (not shown). For spirals with 
large circumference arcs, the low-energy band slopes at the edges of 
the JZ are markedly higher, meaning curvature- dependent Fermi- 
velocities. This dependence is absent in single- 7i band tight-binding 
model since it cannot account for curvature effects. Furthermore, in 
these spirals the chiral operation changes the orbital orientation. The 
angular parts of the p-orbitals can be written as px(l) = px cos(lO), 
Py(l) — Py sm(lO), and pzil) = pz, where / is an integer number for 
consecutive units cells of two non-equivalent carbon atoms. This 
orbital orientation changes the phase of the wavefunction, a 180° turn 
flips its sign; this property modifies the profiles of bonding and anti- 
bonding bands markedly, due to additional cross- coupling terms^^'^^. 
In the simpler spiral example shown on Fig. 3(a), its valence band 
develops minima at X-points and maxima around the F-point. As the 
spiral becomes wider, the electronic structure becomes richer but still 
such peculiar band offset is maintained. This energy uplift and con- 
comitant maximum development around the T point arises because 
the states become more antibonding due to chiral symmetry. 

Zooming over the low energy bands around T point [bottom 
panels of Fig. 3], one can notice that they display intriguing antic- 
rossings close to the Fermi energy as the spirals get wider. This occurs 
due to certain selection rules (similar to those manifested in ordinary 
atomic rings) that the helical states must obey combined with the 
chiral symmetry of the curved space. These anticrossings between 
such chiral branches can be finely tuned by means of external fields. 
This is shown on Fig. 4 which depicts the band structure for the spiral 
ac[l]zz[0-2] while a high intense external electric field is applied 
perpendicularly its axial direction. At such intense electric fields, 
one would expect more impacting changes on the electronic struc- 
ture of the spirals. Nonetheless, such robustness is simply a remark- 
able proof that these helical modes are protected by time- reversal 
symmetry. We can allude the nature of these states by assuming that 
the spirals are composed of q concentric rings embedded in the same 
unit cell, if the spiral respects hexagonal symmetry, i.e. each ring 
contains a number of atoms such as M = 6, 18, 30, 42, . . ., according 
to a tight-binding description, whenever q is even, an energy gap 
opens in the electronic structure (see Fig. 2 in supplementary 
information). Otherwise, metallic edge-states touching the Bril- 
louin zone boundaries will form. The same opening- closing rules 
for the energy gap fail for the first principle results since they take 
into account curvature effects. However the low energy states still 
retain its edge-nature as can be seen from the isosurface plots of the 
local density of states at the Fermi energy for the spiral ac[l]zz[0-2] 
(see Fig. 3 in supplementary material). Ultimately, by comparing the 
results derived from ab initio and tight-binding methods, we suc- 
cessfully confirmed that the electronic structure of graphene spirals 
couples to the helical space backbone. 
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Figure 4 | (red lines) Band structure for the graphene spiral ac[l]zz[0-2] 
(top panel) under the influence of an external electric field of intensity 
0.6 V/A and applied perpendicularly to its axial direction. Horizontal line 
marks the Fermi energy (not shifted to zero). Dotted black lines are the 
energy dispersion for the same spiral for null electric field. 

Discussion 

Previous works have reported that ripples in graphene can be mod- 
eled by coupling the Dirac equation to a curved metric space defined 
phenomenologically from corrugations observed in experiments^^'^^. 
Such covariant formalism gives rise to an effective Hamiltonian 
where it is possible to identify that the electrons on the deformed 
space exhibit a new Fermi velocity. The latter is smaller than the flat 
graphene case. Their findings were confirmed by Raman spectro- 
scopy measurements performed in folded graphene samples^". In our 
spiral examples, electrons propagating along the helical track with 
stronger bending would then manifest smaller effective Fermi velo- 
cities in comparison to those with smoother curvature. This differ- 
ence can be clearly seen on the band structures of 3 where the low 
energy bands associated to each spiral have distinct slopes. Such 
velocity reduction cannot be captured by the single Ti-band tight 
binding approximation since curvature is invisible within this 
description. Another prominent feature that can be distinguished 
from the effective covariant Hamiltonian is the appearance of an 
effective magnetic field pointing perpendicularly to the deformed 
graphene sheet^^'^\ In other words, mechanical strain loaded on 
the distorted section of the graphene mimics the same effects caused 
by an external magnetic field (see supplementary information where 
we demonstrate how the Rashba-like Hamiltonian can be derived 
from an illustrative toy- model example). The field intensity is strictly 
related to how sharp the curvature is. This effective field appears 



naturally from covariant formalism where Dirac equation is solved 
on curved space. Once the metric of the deformed space is known, 
this information is plugged on the Dirac equation and techniques 
resembling to perturbation theory can be used if one assumes that the 
curvature is rather smooth. The solution for this problem can be 
often recognized as the standard graphene Dirac model in the pres- 
ence of an effective potential generated by its own curvature. Such 
method is very efficient in dealing with smooth ripples or light cor- 
rugations on the graphene sheet. Extending such interpretation to 
our systems, one can already expect that the highly intensified curv- 
ature of the graphene spirals will affect enormously their electronic 
response. Since the curvature of the graphene spirals is so remark- 
able, the use of covariant methods following expansion procedures is 
unreliable. In this sense, we must indeed rely on robust ab initio 
methods where all the structures are fully optimized and the effects 
of the curvature are naturally incorporated in the Hamiltonian. 

But even if unreliable, the formalism can help understand the 
origin of the Rashba splitting. In Aharonov-Bohm (AB) devices real 
magnetic fields induce electron wavefunctions an additional phase 
due to the breaking of time-reversal. The phase is proportional to 
magnetic flux penetrating the device geometry and depends on 
whether electron moves clockwise or anticlockwise^^. For the sake 
of simplicity, let us consider a spiral chain with 18 atoms in its unit 
cell (Cig) and investigate its density of states while varying the unit 
cell length or the curvature of its helix [see Fig. 5(b)]. This result was 
obtained within density functional tight-binding (DFTB) formalism 
in which the atomic structure of the system-chain could be fully 
optimized for each unit cell length. Electrons moving clockwise 
(/c+) and anticlockwise (A:—) acquire different phase factors in their 
wavefunctions depending on the local curvature of the spiral. Thus, 
an oscillatory pattern akin to AB interference emerges in the density 
of states [see Fig. 5 (a)] . Only here the AB oscillations arise not from a 
real magnetic field but from a pseudo -magnetic field defined by the 
system topology; the curvature can be seen as a parameter replacing a 
real external magnetic field to produce the same effect. As a result, 
electrons having different pseudospin and being coupled to the 
pseudo -magnetic field can move in opposite directions and still pre- 
serve time- reversal symmetry while exhibiting non-zero AB phase. 

Because of electronic properties similar to topological insulators, 
graphene spirals could naturally be used in quantum computing. 




Figure 5 | (a) Aharonov-Bohm-like oscillations in elongated Cig atomic 
chain being gradually deformed into a spiral shape. Normalized density of 
states vs. energy and elongation of the unit cell (in Angstroms). Line at Ep 
Fermi-level. Horizontal lines highlight the results for some chain 
geometries that are displayed on panel (b). 
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Furthermore, the spiral geometry itself suggests usage as electronics 
components, as nano- solenoids to produce local magnetic fields. We 
can estimate the magnetic field intensity, B, created by electrons 
traveling around the spiral of radius R. By setting the Lorentz force 
and the centripetal force equal, we obtain B = mvp/eR, where m is the 
electron mass, e the electron charge, and Vp the Fermi velocity. 
Assuming R ^ 1 nm and Vj? ~ 10^ m/s (Fermi-velocity of flat gra- 
phene) we obtain an estimate for graphene spiral -generated mag- 
netic field disB = 10^ T. The estimate depends on the Fermi-velocity, 
which our results shows to depend on the precise atomic structure of 
the spiral. The spiral structure can hence be used to sensitively cus- 
tomize the solenoid properties. 

Finally, we discuss the feasibility to realize spirals experimentally. 
We note that natural growth often favors chiral molecules, first and 
most notable example being the DNA. Second example is carbon 
nanotubes grown into helical conformations by using controlled 
creation of pentagon and heptagon defects during the growth pro- 
^ggg53,54 -j-Q synthesize graphene spirals, however, the most relevant 
experiment was the recent controlled synthesis of narrow graphene 
nanoribbons^^. Ribbons were synthesized for controlled ribbon 
widths, controlled edge shapes and —in particular— controlled 
topologies. The control was achieved by organic chemical reactions 
on gold surfaces, where the topology of the nanoribbon was deter- 
mined by the topology of the precursor monomers. Here, we propose 
that graphene spirals could be synthesized by the same controlled 
bottom-up approach. One would only need to choose appropriate 
precursor monomers, preferably from an organic polymer family 
with a helical motif There are also other synthesizing alternatives. 
Namely, the graphene spiral ac[l]zz[0-l] is already a familiar mole- 
cule, the helicene^^. Fairly long helicenes have been synthesized, and 
perhaps related techniques could be extended to synthesize also 
wider graphene spirals^^'^^. Yet another alternative to spiral synthesis 
is to use the viewpoint of array of screw dislocations^^'^^. It was 
reported recently, that in nature closely related structures can be 
found as shown by Rakovan & Jaszczak^^. 

To conclude, the electronic structure of graphene spirals show 
Rashba splitting as a distinct topological signature. The splitting 
can be understood as a consequence of the intrinsic curvature present 
in graphene spirals, as a consequence of the coupling between 
pseudo-spin and the curved helical geometry. The splitting mech- 
anism is similar to the mechanism of band inversion in topological 
insulators^°'^\ The split, low-energy states around F-point are loca- 
lized at the edges and are protected by the spiral topology, being thus 
robust against impurities or lattice distortions. These unique elec- 
tronic properties require neither an external magnetic field nor spin- 
orbit interaction, which is unlike any typical quantum Hall system. 
Therefore, graphene spirals ought to deserve a prominent role as a 
fundamental graphene topology, comparable to the topologies of 
carbon nanotubes and graphene nanoribbons. 

Methods 

Presented band structure calculations were performed using the SIESTA^^ "^" package 
within generalized gradient approximation for the exchange and correlation 
energies^^''^^ Norm-conserving pseudopotentials*^* with relativistic corrections and a 
split-valence double-^ basis of pseudoatomic orbitals with an orbital confining energy 
of 0.05 eV and an energy cutoff of 150 Ry were used, with Perdew-Burke-Ernzerhof 
functionaP^'^^. The k-point sampling contains 6 k-points along the spiral axis 
(simulation cell has length b along spiral axis). Spirals were optimized using 
Hellmann-Feynman forces down to 0.01 eV/A tolerance^^ *^*^. Because of an elastic 
axial stress, the relaxation resulted in an average layer separation of ^ ~ 3.2 A, 
somewhat smaller than graphite interlayer distance. Therefore, because van der 
Waals forces are much weaker than elastic forces, and because their role for the 
electronic structure is insignificant, they could be safely neglected. Also discussion in 
this work weakly might be affected by the fact that super cell optimization shows a 
possible for smallest spiral to change a hexagonal base to a pentagonal. For VASP 
calculations, the projector augmented wave and generalized gradient approximation 
for exchange and correlation energy were used^^ *^^. Kohn-Sham orbitals were 
expanded in plane-wave basis set with energies up to 550 eV and the Brillouin zone 
was sampled over 1X1X4 Monkhorst-Pack grid. For the pseudopotentials. 



Vanderbilt's ultrasoft potentials'^ with cutoff energy of 58 Ry was used. All used 
first-principle methods gave results in excellent agreement. 
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